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ON THE SOUND FIELD OF A POINT-SHAPED SOUND SOURCE 
IN UNIFORM TRAHSLATORY MOTION*^ 

Ry H. Honl 

INTRODUCTION AND SUMMARY^ 


The following report will treat rigorously the excitation of sound 
by a point-shaped sound source (pulsator) in uniform translatory motion 
at subsonic or supersonic velocity through a medium at rest for the two- 
and three-dimensional case (two- and three-dimensional problem) . 

In qualitative respect, the phenomena of sound propagation excited 
by a moving sound source are weXl known. If the construction of surfaces 
of constant phase is based on Huyghens' principle such that the propaga- 
tion in the medium at rest of the elementary waves emanating from the 
sound source is regarded as independent of the momentary state of motion 
of the sound source, characteristic traits' of the sound propagation may 
be understood even on the basis of simple geometric constructions. If, 
for instance, the motion of the sound source visualized as point-shaped 
is uniform translatory along the negative x-axis of a three-dimensional 
Cartesian coordinate system, the sound propagation for a sound source 
moved at subsonic or, respectively, supersonic velocity is represented 
in the known manner by figure 1; according to construction, surfaces of 
constant phase would correspond to the circles plotted. For the motion 
of the source at supersonic velocity, it is, above all, characteristic 
that the excitation of sound always remains limited to the interior of 
a cone - the Mach cone - the (semi) opening angle a of which results 
directly from the construction of the elementary waves (one has 
sin a = c/U when c denotes the sonic velocity, U > c the velocity 
of the sound source) . From the physical point of view, however, there 
arises the question how far the construction of the phase surfaces 
according to the principle of Huyghens* elementary waves can be applied 

*"Uber das Schallfeld einer gleichformig-translatorisch bewegten 
punkformigen Schallquelle . " Annalen der Physik, issue 5, vol. 43, 1943» 
PP* 437-464. 

^Dedicated to Privy Councillor A. Sommerfeld for his 75th anniver- 
sary on December 5, 1943* 

^The present investigation was written within the scope of my work 
for the Kaiser -Wilhelm Institute for Flow Research, Gottingen. 
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rigorously, also the problem of the exact amplitudes or intensities, 
respectively, of the sound field. It will be shown on the basis of the 
rigorous solution of the plane and spatial problem that, on the whole, 
the qualitative description of the sound field according to the prin- 
ciple of elementary waves is justified, but that in particular, especially 
in case of supersonic velocity, characteristic deviations from the ele- 
mentary construction become evident. 

The mathematical theory of sornid propagation will be based on the 
linear propagation equation for the medixmi at rest 



^ = 0 
St2 


( 1 ) 


cp may signify therein for instance the pressure disturbance. ‘It would 
appear rather obvious to utilize the acoustic optical analogies for the 
problem of the sound source in uniform translatory motion and, accord- 
ingly, to apply to (1) a Lorentz transformation where then the sonic 
velocity c instead of the velocity of light would appear in the trans- 
formation formulas as critical velocity. 3 However, this method would 
seem somewhat artificial for the present problem since the transformation 
formulas of the space and time quantities, differently from optics, would 
have no physical significance and their application thus would be a merely 
formal methodical expedient ^ finally, in case of supersonic velocity, one 
would depart even further from the physical starting point, due to the 
different relativity formulations.^ 

We shall therefore adopt below another direct method for solution of 
the problem, making use of the method of the Foiarier integral. Therein 
it is expedient to use a coordinate system with respect to which the sound 
source is at rest at the origin of the coordinate system. Our problem 
then is, except for a Galileo transformation, obviously identical with 
the souind propagation about a sound source in an oncoming flow of sub- 
sonic or supersonic velocity, respectively, assumed at rest at the coordi- 
nate origin. If the oncoming flow has positive x-direction corresponding 


3 

It is known that historically the investigation of the propagation 
equation (l) for the first time gave occasion for setting up linear trans- 
formation formulas for the space and time quantities which are closely 
related to those of the special relativity theory. Compare W. Voigt, 

"liber das Dopplersche Prinzip," Gottingen Nachrichten, I887, p. 4l; see 
also W. Pauli, Enzykl. der Math. Wissensch, vol. 19, p. 543. 

^For the subsonic region, H. G. Kussner, Luftfahrtforsch, vol. I7, 
p. 370, 1940, derived exact formulas for the sound field of a moving sound 
source which agree with the results of this report in an interesting man- 
ner by application of a "Lorentz transformation." 
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to figure 1, the equation of propagation to he solved follows from equa- 
tion ( 1 ) simply hy replacement of the time operator ~ hy — + U ~ 

^t St Sx 

where U is the free-stream velocity. We thus obtain from equation (l) 


If we limit ourselves immediately to the case of a sound excitation 
periodic in time with the frequency co and put accordingly 


cp = u(x,y,z)e"^^ 


the time factor may be split off and we obtain from equation (l') 


Au - -4-1 - ioo + U — \^u = 

<^1 


0 


( 2 ) 


(3) 


Equation ( 3 ) is to be fulfilled in the entire space with exception of 
the coordinate origin at which the sound source is located and of possible 
occurring singular lines or surfaces (delimitation of the Mach cone or, 
for the two-dimensional problem, of the Mach wedge). 

If one sets up an expression for the solution for u in the form of 
a multiple (two or threefold) Fourier integral one has the advantage that 
the coefficients of the representation may be readily given as continuous 
functions in the Fourier space due to the linear character of the differ- 
ential equation ( 3 ); likewise, it is easy to determine simultaneously 
the coefficients of the Foirrier representation in such a manner as corre- 
sponds to the presence of a "point-shaped" sound source at the coordinate 
origin or, for the "plane" problem, of a source distribution visualized 
as rectilinear. It may appear at first as a difficulty of this method 
that the integral representations for the solution u thus obtained due 
to the integrand becoming infinite on characteristic surfaces or lines 
of the real Fourier space (for the spatial or plane problem) are subject 
to indeterminate conditions and without special stipulations regarding 
the com’se of the integration paths in the complex domain have no unique 
meaning. The physical viewpoint, which here as in other propagation 
problems enforces uniqueness, is that aside from the conditions named 
which u must satisfy, a radiation condition must be fulfilled in infinite 
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space. ^ Actually, it is precisely the original indeterminateness of the 
integral which confers the freedom of adjusting the solution of equa- 
tion (3) to the natural physical conditions of the problem (suitably 
manipulating the integration paths in the complex domain) and of deter- 
mining it, in this manner, uniquely. 


1. APPLICATION OF THE METHOD OF THE FOURIER INTEGRAL 


We search for solutions of the propagation equation (5) which cor- 
respond to the standard point source located at the coordinate origin. 
We disregard for the time being fulfillment of a radiation condition. 

We leave the dimension number n indeterminate at first, and put sub- 
sequently n = 2 or n = 3, respectively. 

In order to arrive at solutions of this type in a physically and 
mathematically unobjectionable manner, one will find it expedient to 
start, instead of starting from equation (3), from the inhomogeneous 
differential equation 


Aj^u - icJD + U u = f(xi,X2 • • •) (3a) 

c2 \ axi/ 


in which f(x^^,X2 • • *) signifies a soxarce distribution prescribed at 

first arbitrarily (vanishing sufficiently rapidly at infinity) . Let 
the "unit source" isolated at the zero point now be interpreted as the 
limit of a sequence of continuous source distributions f. Let us there- 
fore first visualize continuous positive functions f which have a large 
maximum at the zero point of the coordinate system (point 0) and are, for 
the rest, to be subject to the condition 



dS = 


1 


ih) 


(dS = dx^ dx^ . . . f is to depend on a parameter in such a 

manner that, to the same extent to which this parameter tends toward a 

^Compare A. Sommerfeld, "Die Greensche Funktion der Schwingungs- 
gleichung." Jahresber. d. DMV. 21, p. 309, 1912. Compare also Frank- 
Mises, "Die Differential- u. Integralgleichungen d. Physik." 
Braunschweig 1955? vol. 2, p. 805. 
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limiting value, f tends toward zero everywhere with exception of the 
point 0. Thus, one has in the lines for the characteristic "delta 
function" 


f = 0 for P 3^ 0 


/f as^ 


1 for P = 0 


(4a) 


where the integral may he restricted to an arhitrarily small neighborhood 
of 0. The formulation of the functions f according to equation (4) 
was chosen for reasons of simplicity. This standardization suggests 
itself physically since there results from equations (3a) and (4) for a 
sound source at rest, thus U = 0, with application of Green's theorem 



(5) 


^ do^ surface element of a sphere of arbitrary smallness surrounding the 

point 0, derivation with respect to the outward directed normal), in 
accordance with the definition of the unit source. 


In order to make the generality of the method for integration of 
equation (3a) in question stand out clearly, one replaces this differ- 
ential equation by 


L[u] = L 






(6) 


where L is assumed to be an arbitrary linear differential expression 
in , . . . with constant coefficients of otherwise arbitrary 

order and dimension nimiber. One now makes for f and u the statement 
of a Fourier integral 


^ ‘Jf ■ ■ da 2 . • . da^^ (7a) 

pn P / \ ^( • • • ) 

”//•••/ da^ • • • da^^ 


u 


(7b) 
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Substitution of equation (jb) in ^ then results 

IT . . . 


= • • 


;i(aiXi+02X2+-.-)dai da2 . . . 


da- 


n 




;i(°^l^l^X2+...)^ dag . . . 


da. 


'n 


and since due to the linearity of L 

a a 


i (aiXi-KLgXg • • •) _ 

L(iai,ia2, . . 


Ij ( ^ — — . , , • \e 


there results from equation (6) by comparison of coefficients 


B^a.| .>8-2^ . . . j - 


A( a^ , ag , 


j^icc-j^, i<X2> • • • ) 


( 8 ) 


(9) 


It now remains to determine /-(a^^ag, • • •) corresponding to the condi- 
tions (1) and (4a). From equation (7a) there follows first by Fourier 
inversion 


A^ ) Eg • • ' ) ~ ^ ^ 


-n 


fe 


“^(^Xp+a2X2-t-. . . )^Q 


and thus in the limiting case (equation (4a)) 


A = (2n) 


-n 


(10) 
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Substitution of equations (9) and (lO) in equation (Tb) finally yields 
the desired solution 

gi(aixi+a2X2+. . . ) 

dttg ♦ • • do^ (11) 

L(iai,ia2, • • .) 

The integral (ll) is identical with the method of N. Zeilon^ for 
preparation of solutions of linear and homogeneous differential equa- 
tions with constant coefficients for a prescribed pole -type singularity 
of the solution. The basic idea of Zeilon's motivation of the integral 
formula (ll) also is the utilization of the identity (8). However, 
whereas Zeilon gives first consideration to the requirement of repre- 
senting the solutions of equation (6) prescribed arbitrary source distri- 
bution of f by means of the "fundamental integral," this integral results 
in our method more directly from Fourier* s integral representation of the 
functions f and u and subsequent limiting process. The method may be 
regarded as mathematically strict if a sequence of suitable analytic func- 
tions f has been selected and the representation of f and u does not 
encounter any difficulties before the application of the limiting process. 
The limiting process to the final formula (ll) may then be performed with- 
out hesitation. 
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We return to the initial equation (l*) and, with the use of L, 
make in it corresponding to equations (2) and (8) the replacements 


St 


— - - > - io> 


— — > la 
8x 


Sy 


-> ip > xy 


Sz 



\ i 

(a2 

f o 

+ p2 
9 

) for 

9 \ 

n = 2 



Li 

+ 

+ 7^j for 

n = 3 

(we write furthermore x, y. 


and 

O', 3 , 7 

instead of 


and . ) Then the fundamental solution u = F becomes according to 
equation (ll): 


N. Zeilon, Arkiv for Matematik 6. 19 H, 9 - 1913 /l^^ compare Frank- 
Mises, Ed. I. S. 862 ff. 

7 nACA editor's note: The original German version of this document 

has this section as number 3> ^ut this is believed to be a typographical 
error and has been changed to number 2 to provide consecutive numbers for 
this translation. 
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(a) For the plane problem (n = 2 ) 


F = - 


i(ox+py) 




da dp 


( 12 ) 


with 


fi = (1 - + p 2 + 2 kMa - 


(12a) 


(b) For the spatial problem (n = 3) 

+ CX3 

1 rrr 


F = 


( 2 n )3 




da dp dy 


(13) 


with 


= (1 - M^ja^ + p2 + ^2 ^ 

if we introduce instead of cd and U the propagation parameter k and 
the dimensionless Mach number M 

k = M = (Ik) 

c c 

We consider first the plane case. Obviously the integrand of the 
fundamental solution (l 2 ) becomes infinite if the Fourier coefficients 
a, P lie on the conic section fig = 0 of the real Fourier plane. The 

equation of the characteristic conic section may be written, according 
to equation (l2a) 

(1 - M2)^/a + — — 'l" + (1 - M^)p2 = k2 (15) 

\ 1 - m2/ 

This equation shows that fig = 0 represents an ellipse (in the limiting 

case M = 0 a circle), parabola, or hyperbola, according to whether 
M<1, =1, or > 1 . Correspondingly, the differential equation (3) is 
of elliptic, parabolic, or hyperbolic type. Figure 2 represents a number 
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of conic sections = 0 for vairious values of the parameter M. All 

conic sections pass through the points a = 0, (3 = -k, and are mirrored 

reflections with respect to the a-aixis. 

For the spatial problem, on the left side of equation (15) is to 

be replaced, in conformity with equation (l5a), by + 7^. The surfaces 
= 0 are generated by rotating the conic sections of figure 2 about 

the ,a-axis and thus represent a sphere (M = O), ellipsoids of revolution 
(0 < M < 1), a paraboloid of revolution (M = l), or bi-sheeted hyperboloids 
of revolution (M > l) . 

Following, both the plane and spatial problem will be treated sepa- 
rately for the elliptic, the parabolic, and the hyperbolic case. The 
elliptic case may be completely traced back to the special case M = 0. 

The parabolic and hyperbolic cases, however, require special considerations 
depending on the various connections regarding the characteristic curves 
^2 = 0 and surfaces = 0, and on the stipulated behavior of the solu- 
tion at infinity. In all cases, the selection of the integration paths 
is the decisive factor. 


3. THE ELLIPTIC CASE (SUBSONIC REGION) 
(a) Plane Problem 


In the two dimensional elliptic case, M < 1, equation (15) represents 
an ellipse. If we write its equation in the normal form 


(a - Oof p2 

a2 1)2 


= 1 


(16) 


we have 


a = 


1 - m2 


b = 


\/l - m2 




kM 

1 - M^ 


(16a) 


If we furthermore make an affine transformation in the a,J3 plane, as well 
as in the x,y plane, by putting 


a' = ^1 - M2(a - Oq) 3' = P I = 


/TT 




Ti = y 


(17) 
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the fundamental Integral (12), (l2a) assumes the form 

+00 


F = - 


,iOo^ 


.i(a'l+p'n) 


(2rt)2 \jl - + p'2 - k2 


da' dp' 


( 18 ) 


where k had been equated 


K 



(17a) 


iOt)X 

Except for the factor ■ , the integral appearing in equation (18) 

VI - m 2 

as a function of the variables I, 'H is therefore the same which would 
result for the case M = 0 with the propagation parameter k = k: . This 
last integral^ however^ represents in the known manner, a cylindric wave 

outgoing from the coordinate origin in since the equation® 


_1_ 

k± 


H^(1 )(kp) 


+00 


i(A|+PT]) 


(at)2JJ ^2 + ^2 


k 2 


dA dP 



(19) 


is valid; therein signifies the first Hankel function of zero 

order if (which will be discussed in more detail later on) the infinite 
integration paths, with respect to A and q provided in equation (l9)^ 
are conducted on suitable paths in the complex A- and ^-plane. (Other- 
wise an incoming or standing cylindric wave would result.) If we regard 
the result of equation (l9) for the present as prescribed, we obtain, if 
we substitute equations (l6a), (l7)^ and (l9) in equation (l8) 


F 


-i^^kx f . \ 

“ti x/i - m 2 “ Vi - m2 / 


®Cf. A. 


Sonmierf eld, elsewhere . 


( 20 ) 
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Turning to the discussion of the terminal formula (20), we remark 

first that the Hankel function for large values of its argument 

kP shows asymptotic behavior 

\ 1 if K p- — ) 

Ho^^^(kp)_> ■' ^ (21) 

iJtKP 

The factor Hq^^^kp) in equation (20) corresponds, therefore, taken by 

itself with consideration of the time dependency selected, compare 
equation (2), to an outgoing cylindric wave in the |, t| plane. The 
space-time dependency of the phase is largely modified by occurrence of 
the factor 


exp 



kx\ 

. r\ I 


The phase ^ thus hecomes, again asymptotically, according to equation (20) 



( 22 ) 


If we consider, for instance, the propagation of the special phase surface 
= - i, an elementary conversion of equation (22) yields 




2 


or, with consideration of equation (l^t-), 

(x - Ut)^ + y^ = c^t^ 


(25) 


This simple result signifies that the sixrfaces of constant phase are propa- 
gated, in relation to the flowing medium, asymptotically at the normal sonic 
velocity c, whereas their center is carried along at the velocity U of 
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the flow, just as to be expected on the basis of Huyghens' construction 
of the elementary waves. Going beyond the asymtotic agreement with the 
geometric construction of the elementary waves, our rigorous formula (20) 
remains correct up to arbitrary proximity to the sound source at any rate 
as long as the idealizations leading to the wave equation (l) ("small" 
amplitudes, "point-shaped" sound source) hold true, thus also in a region 
where Huyghens' principle would fail. 

It is interesting to investigate in addition to the phase the 
behavior of the amplitude. According to equations (20) and (21), the 
decrease of the amplitude takes place asymptotically as l/^. The 
surfaces of constant amplitude are therefore cylindric surfaces of 
elliptic cross-section, the axis of which (parallel z) lies in the sounds 
once visualized as rectilinear. The decrease of the amplitude, and there- 
with also of the intensity, takes place more rapidly in the -x direction 
than in the "ty direction. Here, different from the phase, the direction 
of the flow is never one-sided. 

Finally, the limiting case lines k — ^ 0 will be considered. It 
corresponds to to = 0, thus to a static pressure disturbance, as is 
caused by the flow against a thin rod normally to the flow direction. 
According to equation (l2a) becomes therein 

= (1 - m2)o 2 + p2 (24) 

therewith 

tp i(o3c+Py) 

F i — / / S. da dp 



and after the known transformation 
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Here also the decrease of the pressure disturbance in the tx direction 
takes place more rapidly than in the -ty direction by the factor ^ 


The characteristic appesirance of the argument p =1 


of 




1 - 


^1 - m 2 

p 

+ y instead 


for M = 0 is usually denoted as Prandtl's rule. 


(b) Spatial Problem 


All results derived for the plane elliptic case may be transferred 
directly to the spatial problem. One obtains from equations (l 5 ) and (l 5 a), 
maintaining the substitutions (equation (17)) supplemented by 7' = if) 

^ = z. 


F = - 


,ioox 


- 1-00 


r/T- 


i(a'l+^'Ti+r ) 


da' dS' d7’ 


)/l - a '2 + p '2 + 7'2 . 1,2 


(26) 


(with OQ and k having the same values given in equations (l 6 a) 
and (17a). 


If one, furthermore, makes use of the Fourier representation of the 
spherical wave outgoing from the zero point t = .H = C = 0 


1 

4 jt p 


+00 



—00 


^i(Al+hil+vO 

+ Ii2 + v2 _ 


k 2 


dA dh dv 


P = + T|2 + ^2 


(27) 


(postponing intended remarks on the integration paths), there results 
from equations (26) and (27)^ with consideration of equations (l6a) 
and (18a), the terminal formula 


F = 


-i-iL- 

1 e 

|[l - m 2 


kx 





11 - 


+ y 2 + 



+ y^ + z 2 


(28) 
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The discussion is analogous to that for the plane problem. The 
surfaces of constant phase are spherical surfaces which are propagated 
in relation to the flowing medium at the velocity c, while their center 
is simultaneously carried along. This is rigorously valid up to arbi- 
trary nearness to the sound source. The surfaces of constant amplitude 
and intensity, respectively, are the surfaces p = const; the intensity 

decreases as l/p^, thus more rapidly in the ±x-direction than in the 
ty- and "tz -direction. 


(c) For the sake of completeness, a short supplementary remark 
should be added to the proof of the formulas (19) and (27) for the Fourier 
representation of the cylindrical and spherical wave. Regarding detailed 
treatment of these problems, we refer to A. Sommerfeld, elsewhere. 


Let, for the plane case, the double integral (19) be designated 


+00 




G2 = - 


( 211 ) 2 ':^ 


a2 + ^2 


dX dn 


(19) 


Obviously, this integral has at first not yet a unique meaning, since 

pop 

the integrand becomes infinite on the circle X + d = in the real 
X,U plane. If one introduces in this plane, and likewise in the x,y plane, 
polar coordinates 


X = a cos \lf 


q = a sin 


x2 + n2 = 


I = P cos cp 


T] = P sin cp 




,2 . p2 


and visualizes the integration with respect to the azimuth \lr before 
the integration with respect to p, one may write equation (19)> with 
consideration of the known integral representation of the Bessel function 
of zero order 



r'2jt 


iz cos 



iz 

e 


cos(t-cp) 


(29) 


after a simple transformation, in the form 


G 

2 


Ja(c^p) 


a da 


2it 


0 o2 - 


(19a) 
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In equation (l 9 a) the indeterminateness of the integral is not yet elimi- 
nated due to the integrand's becoming infinite for c = k. The integral 
(19a) represents a cylindrical wave extending to infinity, only when the 
integration path is transferred to the negative -imaginary complex a-plane, 
leaving the point a = k to the left in the manner characterized in fig- 
ure 5* The correctness of the contention is known to result immediately 
from the decomposition of the Bessel function Jq into the two Hankel 

functions and 

Jq(cp) =|^Q(l^(ap) + HQ(2)(ap)| ( 30 ) 

of different asymptotic behavior at infinity of the complex 0-plane, and 
from the deformation of the integration paths indicated in figure 5 . The 
integral (19) may then be reduced to the residue of the constituent part of 
the integrand stemming from the first simamand in equation (50) for a = k 
and yields 


Gg = i Ho^^^(tcp) (19b) 

Because of the asymptotic behavior of large real values of 

the argument (compare equation (21)), equation (l9b) actually corresponds 
to an outgoing (divergent) wave. Had the integration path been trans- 
ferred from equation (l9a) to the positive -imaginary complex o-plane, an 
incoming (convergent) cylindric wave would have resulted, corresponding to 

the asymptotic behavior of . Completely analogous relations prevail 

for the spatial problem (Gj). 

The behavior of the wave outgoing from the zero point divergent 
toward infinity may be comprised according to Sommerfeld into an analyti- 
cal condition, which is denoted as radiation condition. It reads for the 
plane and for the spatial problem, respectively 


lim 

p-)oo 




or, respectively 


lira 



= 0 

(51) 


on the infinitely distant boundary of the plane and spatial region, 
respectively. If G is interpreted as Green's function of the infinite 
domain, the addition of a condition of the type (31) together with the 
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known remaining conditions for G is sirfficient for uniquely fixing 
this function?. 

In contrast to the cylindrical and spherical wave G2 and G^ 

propagated from the zero point isotropically in all directions (compare 
equations (19) and (27)), one deals in the general elliptic case 
according to equations (20) and (28) with unsyrametric radiation. We 
omit formulating for this more general case a radiation condition for 
infinity in analogy to equation (3I) and are content with having reduced 
the general elliptic case 0 < M < 1 to the isotropic problem M = 0. 

The difference between outgoing and incoming wave is therefore brought 
about solely by the selection of the integration path in the general 
elliptic case as well. The possibility of reducing the sound propagation 
in case of subsonic approach flow to isotropic sound propagation is 
physically understandable since the propagation of a sound wave is modi- 
fied by the existence of a flow U < c but, basically, not essentially 
changed; however, the relations become completely different if we now 
t\irn to the hyperbolic and parabolic case. 


4. THE HYPERBOLIC CASE M >1 (SUPERSONIC REGION) 
(a) Plane Problem 


In the hyperbolic case, the characteristic conical section 0^=0 
is a hyperbola. If one writes its equation in the normal form 


(Q^ - Qbf 

a'2 " t'2 


(32) 


one has according to equation (15) 


a' = 


- 1 


b' = 


1/m2 - 


«0 = 


kM 


m2 - 1 


(32a) 


By virtue of the affine transformation 


A 


= /m^ - 1 (a - oq) 


p = S I = 


|/m2 - 


n = y (53) 


9 


Compare also W. Magnus, 


Jahresberichte d. DMV, vol. 52, p. 177^ 
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the fundamental integral (l2) becomes 


F = + 


+00 

1 rr 

^ - l 4 i >,2 - 1*2 . „2 ^ 


( 54 ) 


with 


(3J+a) 

\/m2 - 1 


The position of the characteristic hyperbola 


^2 = -A'" 


+ 


= 0 


in the A,p plane obviously suggests performance of the integration over 
A before that over p (compare fig. 4) . The integration over A then 
leads first to the integral 



the integrand of which has a pole each at A = t|fp^ + k^, in order to 
eliminate the indeterminateness of the integral, the integration path 
leading from -00 to +“ has to be suitably detoured into the complex 
A-plane. The type of transfer of the integration path again can be 
selected only on the basis of physical considerations. Since the soxmd 
wave approaches here supersonic velocity, the sound excitation must be 
required to disappear for x < 0. Since in equation (35) x is replaced 
by I, it is therefore a requirement that the integral (55) yields zero 
for negative values of Except for continuous deformations in which 

no singularities are transversed, the integration path is therewith fixed 
one must make a detour into the negative -imaginary complex A-plane, 
leaving both poles to the left. If | is positive, the integration path 
(compare fig. 5) may be replaced by a path in the infinitely distant 
pos it ive- imaginary A-plane and one circling, in positive sense, each of 

the poles at negative |, in contrast, the shifting to 
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the infinitely distant negative -imaginary Tv-plane results in zero as 
it should. One obtains thus by formation of residues 



. dA = 



for I > 0 
for i < 0 


(35a) 


Therewith the double integral occurring in equation (5^) may be 
reduced, disregarding the factor -2n, for | > 0 to the simple integral 


, +00 


K = 


i^iT) sin 

e dji 




( 36 ) 


2 + k2 


Its calculation is achieved in the following manner. If one puts 


H = K sinh d + y? - k cosh 


dp 


+ ^2 


= dd 


K becomes 


2i > 


''4- CO 


^iK(| coshd+T| sinhd ) -iK(| cosh d-T) sinhd ) 


rdd 


If one substitutes furthermore 


I = T cosh X T) = T sinh X - T]^ = t2 


( 37 ) 


K may, with consideration of the addition theorem 


cosh d cosh X i sinh d sinh X = cosh (d - X) 


also be written 




K = 


iKT cosh (d+x) -iKT cosh (d-x) 


'^1 




kd 


( 36 a) 
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One now uses for the Hankel functions 
representation going back to Heine^*^ 



and 



an integral 





Iz cosh -6 
e 


d-3 





g-iz cosh 






(38) 


Since the appearance of iX beside -d in equation (56a) is obviously 
insignificant, because of the integration limits lying at infinity, one 
obtains therewith 


K 






( 39 ) 


Through the substitution, equation ( 57 ) ^ evidently only a part of 
the entire |,T]-plane is covered; for positive t's, only the region 
characterized in figure 6 by I and cross-hatching, respectively. That 
K shows for negative values of t in the region II, according to 
equation (j6a), also values different from zero, is insignificant, since 
the integral (56) has been defined originally only for positive values 
of It is shown, furthermore, that the integral (56) results in zero 


^*^Heine's integral representations for Hv,h^(z) and 




are 


TT (1)/ \ 2 r” iz coshd a 

H''(z) = (-i) — / e cosh nd di 3 

J Q 


= i^^^ ^ r e"^^°°^^^cosh n^ dd 
^ rt Jq 

with the first representation, for arbitrary n, being valid only for 
the upper (positive-imaginary), the second only for the lower (negative- 
imaginary) z half -plane. Specially for n = 0 both representations are 
valid also on the real z axis, which fact is made use of in the text. 
Compare, for instance, R. Weyrich, "Die Zylinderfunktionen und ihre 
Anwendungen" , B. G. Teubner, 1957 ^ P* 50 * 
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in the regions III and IV. If one puts for this purpose 


I = t' sinh T) = t' cosh V t)^ - ( 373 -) 

(t' > 0 in III, t' < 0 in IV), one obtains from equation (36), with 
consideration of 


sinh ■d cosh \1/ - cosh -d sinh = sinh (d i \(() 


actually 



jgiKT 'sinh(d+it')_gii^''' ’ sinh (d-\|() 


dd 


0 


(59a) 


The function K, defined by the integral (36), thus represents a dis- 
continuous function which assumes values different from zero, only in 
the regions I and II, but disappears in III and IV. 

Finally, one substitutes the expressipn for K in F, equation ( 3 ^)» 
If the factors are combined, there results 


- 1 


(- 2 nK) 


and hence after 


insertion of equations (37) > (59) ^ (52a), and (5^a) 



for |y| ^ ■ — and X > Oj otherwise F = 0 

VM^ - 1 


(J+0) 


The analytical result, that the sound propagation is limited to a 
conical region situated symmetrically to the x-axis, including the positive 
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x-axis^ has the simple physical significance: the appearance of the Mach 

angle a. In fact, there results for the boundary of the propagation 
range for the sound excitation from Huyghens* elementary-wave construction 
immediately (compare fig. 7) 


sin a 


£ = 1 
U M 



X 



(41) 


The wave system represented by equation (40) appears, in contrast to 
the relations in the elliptic case, as a standing-wave system with fixed 
nodal lines ^Bessel function Jq ^ which is modulated by a progressing 

wave (exponential factor) . One can explain this behavior in detail in 
the following manner. If one splits the Bessel function Jq in equa- 
tion (40), corresponding to equation (39), into the two Hankel func- 
tions Hq(^^ and and uses their asymptotic representations^^ 

for suif iciexiLly la.x-ge values of the arguiiient 


KT 



vm2 . 1 


y 


2 


there results 



and thus the space-time dependence of the phase becomes 


(40a) 


4> 


1,2 


_!L_ kx t — 1— ? 2. - 

- 1 - 1 yM^ - 1 ^ 


(42) 


11 


Compare, for instance, 
equation (21a) in the text. 


R. Weyrich, elsewhere, p. 46; compare also 
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with the indices 1 , 2, and the upper and lower sign referring, respec- 
tively, to the origin in the Hankel functions and If 

one now considers for instance the propagation of the phases 0 = + 

an elementary calculation yields just as in the elliptic case 


jt 




(x - Ut)^ + y^ = c^t^ 


m 


The lines (or surfaces) with the constant phase values 


7t 


4 


and 



fill, therefore, at any rate sectors of the same circle 


propagating at the velocity c, the center of which is propagated along 
the positive x-axis at the velocity U > c. It has to be noted that the 
argument kt of the Bessel or Hankel functions must always be positive. 
Hence, there results that, viewed by an observer located at large positive 
X, the first summand (in equation (40a)) represents a convex wave out- 
going at a velocity Wj^ > U, the second summand, a concave wave outgoing 

at a velocity W 2 < U (fig. 8 ). By superposition of the two outgoing 

partial waves, the entire wave system originates. Thus, one may speak 
in the hyperbolic case of directed radiation. 


It should be stressed that the circles moving away, which result 
after construction of the elementary waves, do no longer turn out to be 
curves of constant phase (for rigorous consideration) . Rather, a phase 

shift =-^ exists between the front and rear part of these 

curves (in the sense of the direction of propagation). The gradual phase 
change occurs in the neighborhood of the straight line bounding the Mach 
propagation range where the asymptotic representation of the Hankel 
functions is no longer sufficient (indicated by cross-hatching in fig. 9 ) • 
These deviations from the elementary construction, which only the exact 
theory of the propagation phenomenon can disclose, are comparable to those 
occurring in the theory of the refraction of wave systems. 

According to equation (^0), the amplitude is constant on hyperbola 
branches kt = const, and decreases, according to equation (40a), like 

i/\Jk^ toward the interior of the Mach region. On the straight lines 
bounding the propagation range themselves, there occurs a finite jump of 
F (pressure drop, compression shock) which is the larger the more closely 
U approaches the sonic velocity c. 
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(b) Spatial Problem 

In the transition from the plane to the spatial hyperbolic case 
there results from equations (15) and (l5a) in analogy to equation (34) 
the triple integral 


F = 


1 rfr 

(2rt)^ - 1 


dA djj. dv 


(44) 


with the same significance of Oq^ A^ ... as in equa- 

tions (32a) and (35); furthermore^ we put for reasons of symmetry 7 = 
z = ^. The integration over A corresponds exactly to that of equa- 
tion (35); likewise, the integration path is to be selected in the same 
manner as sub a), with the same motivation as in the case of the plane 
problem. One obtains therewith 





= < 


-2jt 




sin 


for 


for 


^ > 0 
I < 0 


If one now introduces in the 
polar coordinates by putting 

}i = a cos ^ 

T) = p COS Cp 


|i,v plane, as well as in the plane, 

poo 

V = a sin ^ ^ + v = a 

^ = p sin cp t| 2 _|_ ^2 _ p2 


F becomes 


F 


-1 

(2jt)2 |/m2 - 1 



(^lf-9) 


sin ^ Vo^ + 

+ K 2 


a da 


If one now again makes use of the integral representation of the Bessel 
function Jq (equation (29) )^ integrating first with respect to 4^ one 
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obtains further 


F = 


1 f 

- 1 


Jq(po) 


sin i\ja^ + 


a da 


(^ 5 ) 


The remaining^ purely mathematical problem consists in the evaluation 
of the integral in equation (^5)- We shall solve this problem in the 
mathematical appendix (section 6^ integrals with respect to cylinder 
functions) from general viewpoints. Here we should like to remark that 
the integration occurring in equation (45) is easily performed for the 
special cases p = 0 and k = 0 . For p = 0 one obtains, if one sub- 
stitutes \/a^ + = o) the expression which is, at first, indefinite 

Jp-O = r 0 do = f" sin (oil) to (U6) 

Jo 

If one notes that in the integrand of equation (45) the (at first omitted) 
factor Jq(po) has, due to its oscillatory character, a convergence- 

enforcing effect, and that for us only the limiting case of this integral 
for p — >0 is of interest, one obtains therewith (for instance analyti- 
cally, by introducing in the integrand of equation (46a) a convergence - 

enforcing factor e”^^ and finally making the transition to the limits 
p— >-0) in a unique manner. 


^ _ cos 

'Jp=0 f- 


(4?a) 


On the other hand, the integral originating for k = 0 from equa- 
tion (45) 

poo 

'^K=0 J sin (pl) da 

may be easily calculated according to the residuum method, with intro- 
duction of the integral representation of Jq and interchange of the 

integration sequences} the integral represents the discontinuous function^^ 


^^Cf. G. N. Watson, Theory of Bessel Functions, I. Aufl., s. 405, 
equation (6) . 
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Jk=0 



III > p (> 0 ) 
111 < p 


(4Tb) 


Comparison of equations (47a) and (47b) with the corresponding formulas 
for the plane hyperbolic problem suggests generalization to arbitrary p 
and K 


J = 





sin i\ja^ + , 

JJ . — a da 

Va^ + 


f' 


cos 

k |/|2 _ p 2 

i \if 

2 o 2 

V ' 

0 

; - P 


for 

for 


| l | > p (> o ) 

Ul > P 

(48) 


In section 6 it will be proved that the performed generalization is 
actually justified. 

By inserting equation (48) in equation (45) > one may immediately go 
over to the final formula for the spatial hyperbolic casej one obtains 





for y >-=== and x > 0; otherwise F = 0 

I^W2 _ 1 


The discussion of the basic solution (49) in the spatial case is in 
all points corresponding to that for the plane problem. The wave propa- 
gation is restricted to the Mach cone. If one splits the standing wave 

system characterized by the factor cos modulated by the preceding 

complex phase factor again into its complex parts 



the combination of the first summand with the preceding phase factor 
(plane wave) results after the transformation indicated above in an out- 
going convex spherical wave, the combination of the second summand with 
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the preceding phase factor, in contrast, in an outgoing concave spherical 
wave. Altogether there results again directed radiation. Differently 
from the plane problem, however, here the phase on a propagated spherical 
surface 


- U (t - r)j-^ + y2 + z2 = (t - t)^ ( 50 ) 

is exact constant, not only asymptotically constant, since here the phase 
of the two partial systems is exactly defined by decomposition of the 

cos {•••) and obviously no phase shift occurs between the spherical 

zones pertaining to the partial systems. Accordingly, the Huyghens* 
construction of the elementary waves (insofar as it is to contain only 
a statement on the surfaces of constant phase) up to arbitrary proximity 
to the Mach cone is here justified. 



The amplitude is constant on the hyperboloidal surfaces 
x^ 2 2 

— - y'^ - z = const and becomes infinite approaching the Mach 

M^ - 1 

cone. On the x-axis (y = z = 0), the amplitude decreases like l/x the 
intensity thus like l/x^. The occurrence of amplitudes which increase 
with approach to the Mach cone beyond all limits indicates that, rigor- 
ously speaking, the validity range of linear theory has been exceeded. 

In fact, the Mach cone represents a compression shock for a more exact 
representation of which the nonlinear hydrodynamic and thermodynamic 
equations must be used. 


5. THE PAEABOLIC LIMITING CASE: M = 1 

(a) Plane Problem 


The quantity ^2 here becomes according to equation (l2a) 


^2 = 2ka + - k2 (51) 

and after substitution of a = a' + k /2 

~ 2 ko,' + 


( 51 a) 
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The characteristic equation = 0 thus represents in the a',|3 plane 

a parabola which passes through the coordinate origin a' = 0, p = 0 
and is opened toward the negative a' -axis. The fundamental solution F, 
equation (l2), may be written 


F = 


1 ^2^ 
e 


(2n)‘ 


-f-00 


dp e 


ipy 



ia'x 

e 

2ka' + p2 


(52) 


The integration over a' may now easily be carried out according to the 
residuum method. Here again the choice of the integration path is decisive. 
We select it in such a manner that we integrate in the complex a '-plane from 

-00 to +00 on leaving the pole of the integrand at ex' =- ^ to the left. 

2k 

Then we may replace the integration path, for positive values of x, by a 
circling of the pole a' = in positive sense whereas for negative 

values of x, there results zero. Thus 




da' = s 


-« 2ka' + 


2j(i 

2k 

0 


-iJ^; 


2k 


for X > 0 
for X < 0 


(55) 


This result justifies subsequently our choice of the integration path. 

In fact^ it is to be expected according to the principle of the elementary 
waves that for a sound source approached by a flow at the velocity U = c 
no sound excitation is brought about for x < 0 since the sound propaga- 
tion is carried along by the flow at the velocity c in the direction of 
the positive x-axis. 


Substitution of equation (55) into equation (52) now yields for 

X > 0 


F = 


jti 

(2rt)2k 


i|x 

e 



(52a) 


By quadratic supplement of the exponent, one readily finds hence 

j_k x^+y^ 




, +00 ^ P 

e da 


(52b) 
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The integral appearing here is well known as "Fresnel's Integral" 



Therewith, one finally obtains 


r 


e'^ e 2 X 


F =^2>/ar 

0 




for X > 0 
for X < 0 


( 5 ^) 


First, one considers again the curves (or siorfaces, respectively) of 
constant phase. If one inserts in the expression for the phase 


il> + = $ ' = yT 

4 

there results because of k = u/c immediately 

l^x - c (t - t) ^ ^ + y^ = c^ (t - t)^ ( 55 ) 

Thus, the surfaces $ = const (t = const) form, as is to be expected 
according to the construction of the elementary waves, a system of circles 
which are tangent at the coordinate origin x = y = 0 and the centers of 
which fill the positive x-axis (fig. 10), The propagation of the phase 
surfaces in the mediimi takes place to all sides and up to arbitrary proxi- 
mity to the sound source at the velocity c. On the other hand, it is 
noteworthy that the decrease of the amplitude over the entire front occurs 

independently of y with l/ \Tx. With the approach of positive x-values 
to X = 0 the amplitude increases beyond all limits; for x < 0 it is 
zero. 
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(b) The Spatial Problem 

becomes according to equation (ija-) and with replacement of a 
by a* as in the plane case 


ft, = 2ka + = 2ka* + p^ + 7^ 


(56) 


therewith the fundamental integral (13), if one visualizes the integration 
over a' carried out in a manner exactly corresponding to equations (52) 
and (55) 


F = - 


Xjx nr 


(2n)5 


i(Py+7z) p +00 
dp d 7 e / da* 


la'x 


2ka* + P^ + 7^ 


i^py+7Z- 


^ ^ ^ .0. .>0 

(2rt)5k 

for X < 0 


(57) 


If one uses furthermore polar coordinates in the p,7 plane just as in 
the y,z plane 


P = 0 cos 
y = p cos Cp 


OOP 

7 = 0 sin ^ P^ + 7 ^^ = 0 


z = p sin cp 


y^ + = p^ 


(58) 


one may write the integral (57) fox' x > 0 

00 2 

// 


■P ’xx 
F e 


V 00 2 Jt y V 0^ 

i|x p n iop cos (i-cp)-!— x 


(2n)5k 


o d 0 d\|r 


00 


(57a) 


First performing the integration over \|/, one obtains with consideration 
of the integral representation of the Bessel function Jq equation (29) ^ 

from equation (5?a) 
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4rtk 

The integral which appears here is 

H. Weber^2 



jQ(pa)e a da (5Tb) 

of the type of the integral named after 



-h^t^ 

e t dt 



e 


a^ 

l+b^ 


and may therefore be reduced to an elementary function; for a = 

15 

1 


b2 = 


2k 


there results 



(po-) 


-i^2 


0 da 


ix 


ik £2 
2 X 


and finally 


F 



,-k x^+P^ 
2 ^ 


e 

X 


for 


X > 0 


0 


for X < 0 


(59) 


The structure of this formula is perfectly analogous to that of 
equation ( 5 ^) for the plane problem. The surfaces of constant phase are 
spherical surfaces which are all tangent at the point x = p = 0 and the 
centers of which again fill the positive x-axis. The amplitude is con- 
stant on the planes x = const (thus independent of p) and, corresponding 
to the intensity relations for the spatial problem, decreases for positive 
X as l/x. 

Furthermore, the analogy between equations (54) and (59) and the 
basic solution (l9b) of the plane problem (compare also the asymptotic 
representation (equation (21a)) and of the corresponding spatial problem 


N. Watson, Theory of Bessel-F-unctions, Cambridge 1922, p. 592. 
^^Convergence of the integral exists as long as |arg b| = ^* 
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for the case M = 0 of disappearing approach flow is remarkable. Finally ^ 
it should be noted that the formulas (5^) and (58) as limiting cases 
M — » 1 can be represented only from the hyperbolic case M > 1, not the 
elliptic M < 1. 


6 . MATHEMATICAL APPENDIX 
Integrals Over Cylinder Functions 


It remains to append the proof for the integral formula (A8) . Since 
this proof requires a few more general considerations, we shall, in con- 
nection with it, derive a few more related integral formulas. 

We start from Heine-Schafheitlin' s integral representation ( 38 ) of 
the Hankel functions and which we write 


■'+<» iz cosh 


d-a = iitHo^^^(z) 


+00 -iz cosh *5 


da = -iitHo^^^(z) 


(58) 


We now consider the integrals 




^-00 +1^1^ +K I +ijiT| 






d|i Ko = 


J^CO 


+iiaq 
e dp 


( 60 a, b) 


As on page I 8 we substitute 

yj[x^ + = K cosh a p = K sinh a 


dp 




= da 


I = T cosh X T1 = T sinh X ^ = 


2 ^ ^2 


(57) 


and thus obtain, according to the addition theorem of the hyperbolic 
functions and to the integral formulas ( 38 ) 


- f 


+00 +1kt cosh (a+X) 


Ki = / e 

^-00 


da = i^Hq 


( 1 ) 


(k _ t)2^ 


(6la) 
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> +00 


Kg - 


^-iKT cosh(^-X) ( 3^3 - T)^) 


(6lb) 


both formulas are valid for the range of representation of the coordi- 
nates Tj by the equation (57) real), that is, for \i \ ^ |ti| . 

In exactly the same manner as on page 20^ one proves that' Kp 
and K 2 disappear for ||| < It |* integrals ( 60 a) and ( 60 b) thus 

represent discontinuous functions which assume values different from 
zero only in the regions I and II of figure 8, but disappear in the 
regions III and IV. Thus we have 




_ Tl^) 


for 

for 


III ^ hi 

III < 1 T |1 


(62a) 



-i 

e 







. - n^) for 

0 for 


Ul > Ini 
III < Ini 


(62b) 

In order to indicate the discontinuous character of the integrals 

and K 2 in the notation, we shall write below instead of the right sides 

of equations ( 62 a) and ( 62 b) abbreviatedly iitHQ^^^ j 

and ) , indicating by a bar above the symbol that the 

Hankel functions are to be applied only for |^| ^ however, zero. 

We now use a new integration variable v in equations ( 62 a) and ( 62 b) 
by putting 
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We thus ohtain 


P +00 o 0 p + 

— / • • • “ / • • • ^ I • • • ^ [ 1 


0 


+00 )fv^ - 1^2 


iv^ 


dv + / S 

K 


+“ _+i\/v^-K^n 


/v^ - 


dv 


or else 

Ki = 2 


^ |/v2 . k2 


+iv^ , . „ 

e dv = irtH 


- 1 ^) 


(63a) 


In exact analogy 


K2 = 2 J dv . -1 «Ho< 2)(, ipT^) 

FV^ - 


(65b) 


By addition and subtraction of equations (63a) and (63b) there results 
with consideration of the mutual connection of the cylinder functions 

r 't r 

Jq(^) = I + Ho^^^z) “ No(z) = Ho(l)(z) -Hq(2)(z) 

^Jq Bessel function, Nq Neumann function of zero order then further 


|\/v^ - 

/v2 - k .2 


cos v| dv = 


^ ” cos r[]lv^ - j n ~ ( \U^ 2^1 

JJ . sin v| dv = + I 

K VvS _ k2 


( 64 b) 
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(the bar over the symbol having the same significance as above). By 
application of Fourier's inversion formulas 


g-i (u) = f fq(t) cos ut dt 

-^0 

fq(t) 

_ 2 

P 00 

/ Si(u) cos ut du 

’ 0 

r “ 



p 00 

gp(u) = / fp(t) sin ut dt 


2 

TT 

/ g 2 (u) sin ut du 

J 0 


Jl 

' 0 


one now obtains from equations ( 64 a) and ( 64 b) 




If we finally put 



the integral formulas (65a) and (65b) assume the form 
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The first of the integral formxolas (66) is identical with equation (48) 
as one recognize'fe immediately if one makes in the designation the 
replacements 


K ^ p V > i T > a T] — > K 

quod erat demonstrandum. 


Translation by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 




moved at subsonic velocity and at supersonic velocity. 
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Outgoing convex wove 
Outgoing concave wave 



Figure 8.- Splitting of the wave system into two partial systems. 





